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INTRODUCTION 


Surface and corner cracks can occur in many structural components. These 
cracks can cause premature failure of landing gear of aircraft, spars, 
stiffeners, and other reinforcements. Accurate stress analyses of these crack 
components are needed for reliable prediction of crack-growth rates and fracture 
strengths. 

Surface and corner cracks in plates have received considerable attention in 
the literature in the past sixteen years [1-20]. Various methods: finite- 
element method with singularity elements [3,4,8,10,12,13,17,19], finite-element 
method with displacement hybrid elements [5,6], alternating method [1,2, 7, 9], 
finite-element-alternating method [14,15,18,20], and boundary-integral equation 
method [11], were used to obtain stress-intensity factors for these cracked 
components. Stress-intensity factor equations were obtained by fitting 
empirical equations to the stress-intensity factors obtained from numerical 
analyses [16,17]. 

Surface and corner cracks at holes have received much less attention than 
surface cracks in plates. Smith and Kullgren [9], and Raju and Newman [12], 
analyzed corner cracks from a circular hole by the alternating method and the 
finite-element method (with singularity elements), respectively. Nishioka and 
Atluri [15,18] analyzed corner cracks from a circular hole and corner cracks in 
a lug by the finite-element-alternating method. The stress-intensity factors by 
all of the above mentioned methods agreed well with one another [12,15], except 
in the region where the crack intersected the hole boundary. In this region, 
the stress-intensity factors from Raju and Newman [12,16] showed a precipitous 
"drop-off". The stress-intensity factors calculated by other investigators (see 
refs. 12 and 15) did not show the large drop-off near the hole-crack 
intersection. Thus, the primary purpose of this paper is to investigate the 
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reasons for the large drop in the stress-intensity factors where a surface crack 
meets a hole or notch boundary. In this investigation, an evaluation of the 
previously developed finite-element models used by Raju and Newman [12,16] and 
of the methods used to determine the stress-intensity factors was performed. 

The second purpose of this paper is to present stress-intensity factors for 
a semi -ell iptical surface crack located at the center of a semi-circular edge 
notch in a finite thickness plate. This crack configuration has not been 
analyzed in the literature. The plates were subjected to remote tension loads 
and a wide range in configuration parameters were considered. The ratio of 
crack depth along the notch root to crack length away from the notch root (a/c) 
ranged from 0.4 to 2; the ratio of crack depth to plate thickness (a/t) ranged 
from 0.2 to 0.8; and the ratio of notch radius to plate thickness (r/t) ranged 
from 1 to 3. The r/t ratio of 3 was chosen because this particular 
configuration was tested in an AGARD study on short-crack growth behavior [21]. 


SYMBOLS 

a half-depth of surface crack along notch root 

b width of plate 

c length of surface crack from notch root 

F stress-intensity boundary-correction factor 

sn 

h half-height of plate 

K stress-intensity factor (mode I) 

Q shape factor for elliptical crack 

r radius of the semi-circular notch or circular hole 

S remote uniform tensile stress 

t half-thickness of plate 
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x,y,z Cartesian coordinate system 

v Poisson's ratio 

parametric angle of ellipse 

THREE-DIMENSIONAL FINITE-ELEMENT ANALYSIS 

A three-dimensional finite-element analysis was used to calculate the mode 
I stress-intensity factor variations along the crack front for a surface crack 
emanating from the center of a semi-circular edge notch (or a circular hole) in 
a plate subjected to tensile loading, as shown in Figure 1. In these analyses, 
Poisson's ratio was assumed to be 0.3. 

Due to symmetry, only one quarter of the specimen was modeled for the 
notch configuration. Symmetry conditions on the x = 0 plane were used to model 
the circular hole configuration. A typical model is shown in Figure 2. This 
model has 8 wedges on the crack plane and 4 layers on the hole. Six-noded 
pentahedral singularity elements were used along the crack front and eight-noded 
hexahedral elements were used elsewhere. The finite-element models were 
subjected to remote tensile loading on the y = h plane. Stress- intensity 
factors were calculated using a nodal-force method. Details of the formulation 
of these types of elements and the nodal -force method are given in References 8, 
10 and 12. Details on the development of the finite-element models will be 
discussed later. 

Stress- Intensity Factor 

The remote tensile loads cause only mode I deformations. The mode I 
stress-intensity factor K for any location clong the crack front was taken to 
be 


3 



K = S (wa/Q) 


( 1 ) 


F sn < a//t ’ a/c ’ r / t ’ * ) 

The half-height of the plate, h, and width, b, were chosen large enough to have 
a negligible effect on stress-intensity factors (h/b > 2 and r/b =* 0.05). 
Values of F , the boundary-correction factor, were calculated along the crack 
front for various combinations of parameters (a/t, a/c, r/t and <f> ). The crack 
dimensions and parametric angle, <f> , are defined in Figure 1. The shape factor 
for an ellipse, Q, is given by the square of the complete elliptic integral of 
the second kind. Empirical expressions for Q (taken from ref. 10) were 


Q = 

1 + 1.464 (a/c) 1 ' 65 

for a/c < 1 

(2a) 

Q - 

1 + 1.464 (c/a) 1,65 

for a/c > 1 

(2b) 


Evaluation of Finite-Element Models and Methods 
Stress-intensity factor distributions in References 12 and 16 for surface 
and corner cracks at a hole showed a precipitous "drop-off" near the region 
where the crack meets the hole boundary. This drop-off was believed to have 
been caused by the "boundary-layer" effect. However, comparisons (see refs. 12 
and 15) with other solutions in the literature showed that the stress-intensity 
factors calculated by other investigators did not exhibit this large drop-off 
near the boundary-layer region. Therefore, a study was carried out to determine 
whether this drop-off is an artifact of the modeling or is a real behavior. The 
models used in References 12 and 16 were re-examined. This examination revealed 
that the procedures used to generate the solids produced models with some "ill- 
shaped" elements in the region where the crack intersects the hole boundary. 
These ill-shaped elements were produced when the "wedge" region along the crack 
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plane was connected directly to the "layer" region to model the hole or notch. 
For r/t - 1, the elements in the region where the crack meets the hole boundary 
had aspect ratios (ratio of largest to smallest dimension) of about 150. The 
aspect ratio of these elements increased linearly as the r/t ratio increased, 
that is, for r/t = 3 the aspect ratio was 450! 

The "ill -shaped" elements were defined as those elements that had aspect 
ratios greater than 100. These elements do not have any influence on the 
results if they are in regions of small stress gradients. However, they tend to 
produce a stiffer model if they are located in regions of large stress 
gradients. Because these elements were found in the region where the crack 
intersects the hole boundary, an investigation was undertaken to study the 
influence of these ill -shaped elements on the stress-intensity factor solutions. 
Three different finite-element models were developed and two different methods 
were used to extract the stress-intensity factors from the finite-element 
solutions. 

Models. - Model A used finite-element patterns that were similar to those 
used in References 12 and 16 with the "ill -shaped" elements near the hole-crack 
intersection. Model B was a newly-developed finite-element model similar to 
Model A but without ill -shaped elements. Both Models A and B used a polar mesh 
and singularity elements around the crack front, as shown in Figure 2. Model C, 
on the other hand, used a rectangular mesh and non-singular elements around the 
crack front. Also, Model C did not have any ill -shaped elements near the hole- 
crack intersection. In Models B and C, the ill-shaped elements were avoided by 
not directly connecting the wedge region to the layer region near the hole-crack 
intersection. A small offset (less than 0.04r) was used before the layers were 
added to model the hole or notch configuration (see fig. 2). The x-coordinates 
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of the nodes in the "offset" region were adjusted so that they lie on the 
circular boundary. 

Two different in-house (NASA Langley) computer programs, discussed in 
References 12 and 23, were used to analyze the various models. The different 
computer programs were used to verify the consistency of the results and to 
determine whether any undetected errors existed in the computer codes and 
analyses. 

In Models A and B, pentahedron singularity elements were used all along the 
crack front and 8-noded hexahedron elements were used to model the rest of the 
solid. Furthermore, a 3x3x3 Gaussian integration scheme was used to form the 
element stiffness matrices. In Model C, non-singular 8-noded hexahedron 
elements were used everywhere to model the solid. A 2x2x2 Gaussian integration 
scheme and reduced integration for the shear terms [23] were used in conjunction 
with Model C. In contrast to Models A and B, where the models had a polar-mesh 
pattern around the crack front, Model C had a ^ectangul ar-mesh pattern around 
the crack front. 

Methods . - Two methods for extracting stress-intensity factors from the 
solutions were also studied. The methods were the nodal-force method [8,10,12] 
and the virtual -crack-closure technique (VCCT) [23,24]. In the nodal -force 
method, the forces normal to the crack plane and ahead of the crack front were 
used to calculate the stress-intensity factors. This method does not require 
the assumption of either plane stress or plane strain to obtain the stress- 
intensity factor. The VCCT calculates the strain-energy release rate using the 
forces normal to the crack plane and ahead of the crack front and the relative 
displacements of the crack faces behind the crack front. The stress-intensity 
factor is then obtained from the calculated strain-energy release rate by 
assuming either plane-strain or plane-stress conditions. Plane-strain 
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conditions were used everywhere along the crack front except where the crack 
meets the notch surface {<j> = 90°). At this location, plane-stress conditions 
were used. These types of assumptions along the crack front are widely used in 
the literature [4,11,13,23,24]. 

Comparison of stress-intensity factors. - Three specimen configurations with 
r/t = 1 or 3 were chosen in the evaluation study. All configurations had a 
semi-circular (a/c = 1) crack with an a/t ratio of 0.2. One configuration had a 
surface crack emanating from a circular hole (r/t = 1) while the other 
configurations had a surface crack emanating from a semi-circular edge notch 
(r/t = 1 and 3). All models had 14 wedges to model the crack plane and 10 
layers to model the notch. Models A and B had about 15,000 degrees-of-freedom 
(dof) while Model C had about 18,000 dof. 

Figure 3 presents the stress-intensity factor distributions along the crack 
front from Models A, B and C for a surface crack located at the center of a 
circular hole with r/t = 1 . The results from Models B and C agreed well (within 
1.5 percent) with each other all along the crack front. Model A, with the "ill- 
shaped" elements (circular symbols [16]), gave stress-intensity factors that 
were lower than the results from the other models. The discrepancy was quite 
large for <j> > 60°. The results from Models B and C do not show the large drop- 
off exhibited by Model A. 

The solid curve is an empirical equation for a surface crack at a circular 
hole that was based on previous finite-element results [16]. Comparison with 
the present finite-element results suggests that the equation is still accurate 
(within 5 percent) for this crack configuration although a slight modification 
would be necessary near the hole-crack intersection. 

Figures 4 and 5 present stress-intensity factors for the same crack shape 
and size as shown in Figure 3 (a/c=l; a/t=0.2) but for a surface crack at a 
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semi-circular edge notch with r/t ratios of 1 and 3, respectively. The nodal- 
force results from models with and without ill -shaped elements (Models A and B, 
respectively) differed all along the crack front especially for the results with 
r/t = 3. Again, the largest discrepancies occurred for 4 > > 60°. However, the 
nodal -force results (Model B) were in good agreement with results from Model C 
using the VCCT. The maximum difference between these two sets of results was 
only about 2 percent. 

Again, the solid curves in these figures show stress-intensity factors 
calculated from an empirical equation [25]. This equation was previously fitted 
to finite-element results for a surface crack at a hole [16] but adjusted for 
the stress concentration differences between a hole and a notch. Comparison 
with the present finite-element results suggests that the equation is still 
accurate for this crack configuration although a slight modification would, 
again, be necessary near the hole-crack intersection (maximum difference at this 
location was about 5 percent). However, for r/t = 3, the equation consistently 
underestimated the stress-intensity factors (5 to 10 percent). A reason for the 
lower estimates is probably because the equation was extrapolated to r/t values 
beyond those used in fitting the equation (0.5 < r/t < 2). Thus, a modification 
to the equations is required for crack configurations with the r/t ratios 
greater than 2. 

Discussion. - Recall that the nodal -force method was used to extract the 
stress-intensity factors from the models with and without the ill-shaped 
elements (Model A and B, respectively). On the other hand, the VCCT was used in 
conjunction with Model C. All models without the ill-shaped elements gave 
nearly the same stress-intensity factors all along the crack front. This 
suggests that the two methods of extracting stress-intensity factors are 
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reliable, provided that the models are well behaved. The models with the "ill- 
shaped" elements tended to make the model too stiff and, hence, gave lower 
stress-intensity factors than models without these elements. Therefore, 
improved finite-element models that were similar to those used in References 12 
and 16 but did not contain any ill -shaped elements were developed for a wide 
range in configuration parameters. These new models, in conjunction with the 
nodal-force method, were used to analyze a surface crack at a semi-circular edge 
notch herein. These models had about 15,000 dof. The results obtained with 
these models are discussed in the following section. 


STRESS- INTENSITY FACTOR SOLUTIONS FOR SURFACE CRACK AT SEMI-CIRCULAR EDGE NOTCH 
A surface crack emanating from a semi-circular edge notch in a plate 
subjected to remote tensile loading was considered herein. A wide range in 
crack sizes (a/t), crack shapes (a/c) and notch sizes (r/t) were analyzed. The 
ranges in crack and notch parameters were: 

a/t = 0.2, 0.5, and 0.8 
a/c = 0.4, 1, and 2 
r/t = 1, 2, and 3 

These particular crack configurations were chosen to cover the range of crack 
shapes and sizes that have been observed to grow in experiments. Note that the 
r/t = 3 configuration corresponds to a specimen used in an AGARD study on short 
cracks (see refs. 21 or 25). For each combination of parameters, stress- 
intensity factors were obtained all along the crack front. 


9 



The normalized stress- intensity factors obtained with the three- 
dimensional finite-element analyses are presented in Tables 1-3. Typical 
results are shown and discussed herein. 

The normalized stress- intensity factor variations all along the crack front 
for the AGARD short-crack specimen with r/t = 3 are shown in Figures 6 and 7 for 
various crack shapes (a/c = 0.4 and 1, respectively). In these figures, the 
inserts show part of the crack plane with the relative size of the notch, half- 
thickness of the plate, and the various crack shapes and sizes considered. For 
both crack shapes, the maximum normalized stress-intensity factors occurred near 
the intersection of the crack front and the notch boundary ( <j> = 90°). This 
trend is expected because the crack front lies in a region that is influenced by 
the notch. Hence, surface cracks would be expected to grow more along the bore 
of the notch than away from the notch and approach an a/c ratio greater than 
unity. Experiments on surface cracks at notches [21,25] tend to support this 
observation. 

Figure 8 shows the normalized stress-intensity factor variations along the 
crack front for a semi -ell iptical surface crack with an a/c ratio of 2 and r/t = 
1. In this case, the normalized stress-intensity factors are largest at the 
maximum crack length from the notch root (tf> = 0°). For deep cracks (a/t = 0.8), 
the normalized stress- intensi ty factors are nearly constant all along the crack 
front . 

The variation of normalized stress- intensity factors all along the crack 
front for various r/t ratios with a/c = 0.4 and a/t = 0.8 is shown in Figure 9. 

As expected, larger values of notch radii (hither r/t ratios) gave higher values 
of normalized stress-intensity factors for the same crack shape and size. This 
trend is expected because the crack surfaces c're subjected to a higher normal 
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stress component for the larger r/t values. Similar trends were observed for 
all crack shapes and sizes analyzed and, hence, the results are not shown. 

Figures 10 and 11 present the normalized stress-intensity factors at two 
locations along the crack front as a function of a/t for two different crack 
configurations. The stress-intensity factors at <t> = 0° and near the end of 
the other axis ( <f> = 80°) are presented. The 80° location was selected to avoid 
the boundary-layer region. Finite-element results were obtained for a/t ratios 
of 0.2, 0.5 and 0.8. The limiting solutions for a surface crack at a notch as 
a/t approaches zero were estimated from stress-intensity factor equations for a 
surface crack in a plate [16] and a stress concentration factor of 3.17 [25]. 

The normalized stress- intensity factors at <f> = 0° and 80° are shown in 
Figure 12 as a function of notch size (r/t) for an a/c ratio of 0.4. The 
results for the shallow crack (a/t = 0.2) are shown as solid curves, whereas the 
results for a deep crack (a/t = 0.8) are shown as dashed curves. These figures 
show that larger size notches (higher r/t ratios) tend to give higher normalized 
stress-intensity factors at both ends of the crack. 


CONCLUDING REMARKS 

A re-evaluation of the finite-element models and methods used to extract 
stress-intensity factors from finite-element solutions for surface cracks at 
stress concentrations was performed. Previous finite-element models used by 
Raju and Newman for surface and corner cracks at holes were shown to have "ill- 
shaped" elements at the intersection of the hole and crack boundaries. These 
ill-shaped elements tended to make the models too stiff and, hence, gave lower 
stress-intensity factors near the hole-crack intersection than models without 
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these elements. The improved models employed 8-noded elements and stress- 
intensity factors were calculated by both the nodal -force and virtual -crack- 
closure methods. Both methods and different models, without the ill -shaped 
elements, gave essentially the same results. 

Comparisons made between the previously developed stress-intensity factor 
equations and the results from the improved models agreed well (within 5 
percent) for configurations with notch-radii-to-plate-thickness (r/t) ratios 
less than 2. Larger discrepancies were shown for an r/t ratio of 3. Here, the 
stress-intensity factor equations gave results that were 5 to 10 percent lower 
than the results from the improved models for a semi-elliptical surface crack at 
a semi-circular edge notch. 

Stress-intensity factors for semi -ell iptical surface cracks emanating from 
from the center of a semi-circular edge notch in a plate subjected to remote 
tensile loading were obtained by three-dimensional finite-element analyses. A 
wide range of crack shapes, crack sizes, and notch-radius-to-plate-thickness 
ratios were considered. The range of crack sizes, defined by crack-depth-to- 
plate-thickness (a/t), considered was 0.2 to 0.8. The range of crack shapes, 
defined by the ratio of crack length along the notch to the crack length away 
from the notch root (a/c), considered was 0.4 to 2. The range of notch sizes, 
defined by the notch-radius-to-plate-thickness (r/t), considered was 1 to 3. 

The configuration with r/t = 3 was chosen because this is the configuration that 
was tested in an AGARD cooperative test program on growth of short cracks at 
notches . 
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Table 1. Boundary-correction factors, F $n , for semi -elliptical surface 
crack emanating from the center of a semi-circular notch in a 
plate under tension. (r/t-1; F $n * K/[S(>ra/Q)^]) 


a/c 

^ (deg ) 


a/t 


0.2 

0.5 

078 

0.4 

0.0 

1.113 

0.869 

0.833 


22.5 

1.284 

0.997 

0.967 


45.0 

1.678 

1.296 

1.274 


67.5 

2.190 

1.733 

1.774 


78.8 

2.567 

2.129 

2.232 


84.9 

2.809 

2.445 

2.631 


88.0 

2.938 

2.655 

2.951 


90.0 

2.727 

2.548 

2.944 

1 

0.0 

2.402 

1.865 

1.697 


22.5 

2.442 

1.904 

1.716 


45.0 

2.574 

2.046 

1.841 


67.5 

2.826 

2.355 

2.211 


78.8 

3.043 

2.646 

2.599 


84.9 

3.192 

2.858 

2.898 


88.0 

3.263 

2.975 

3.089 


90.0 

2.858 

2.656 

2.826 

2 

0.0 

1.988 

1.545 

1.299 


22.5 

1.942 

1.513 

1.245 


45.0 

1.831 

1.453 

1.173 


67.5 

1.741 

1.456 

1.294 


78.8 

1.719 

1.491 

1.409 


84.9 

1.748 

1.546 

1.504 


88.0 

1.781 

1.590 

1.566 


90.0 

1.632 

1.467 

1.431 
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Table 2. Boundary-correction factors, F $n , for semi-elliptical surface 
crack emanating from the center of a semi-circular notch in a 
plate under tension. (r/t=2; F $n - K/[S(wa/Q) ) 


a/c 

(deg) 


a/T 


0.2 

075 

O 

0.4 

0.0 

1.405 

1.081 

0.989 


22.5 

1.613 

1.245 

1.156 


45.0 

2.070 

1.634 

1.553 


67.5 

2.596 

2.174 

2.207 


78.8 

2.943 

2.595 

2.764 


84.9 

3.151 

2.880 

3.182 


88.0 

3.252 

3.044 

3.469 


90.0 

2.982 

2.844 

3.347 

1 

0.0 

2.789 

2.316 

2.128 


22.5 

2.822 

2.357 

2.152 


45.0 

2.929 

2.497 

2.303 


67.5 

3.138 

2.780 

2.715 


78.8 

3.323 

3.028 

3.105 


84.9 

3.454 

3.201 

3.378 


88.0 

3.509 

3.285 

3.533 


90.0 

3.052 

2.885 

3.161 

2 

0.0 

2.309 

1.832 

1.595 


22.5 

2.247 

1.786 

1.523 


45.0 

2.098 

1.683 

1.413 


67.5 

1.975 

1.627 

1.506 


78.8 

1.945 

1.624 

1.594 


84.9 

1.978 

1.662 

1.670 


88.0 

2.015 

1.697 

1.720 


90.0 

1.847 

1.554 

1.554 
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Table 3. Boundary-correction factors, F sn , for semi-elliptical surface 
crack emanating from the center of a semi-circular notch in a 
plate under tension. ( r/t =3 ; 


sn 


= K/[SUa/Q) 1/2 ]) 


a/c 

*(deg) 


a/t 


0.2 

0.5 

0.8 

0.4 

0.0 

1.603 

1.253 

1.131 


22.5 

1.831 

1.442 

1.325 


45.0 

2.318 

1.882 

1.788 


67.5 

2.839 

2.460 

2.534 


78.8 

3.167 

2.878 

3.135 


84.9 

3.361 

3.145 

3.559 


88.0 

3.452 

3.290 

3.834 


90.0 

3.151 

3.044 

3.655 

1 

0.0 

3.017 

2.592 

2.429 


22.5 

3.045 

2.630 

2.454 


45.0 

3.138 

2.761 

2.610 


67.5 

3.328 

3.022 

3.032 


78.8 

3.503 

3.251 

3.421 


84.9 

3.629 

3.411 

3.685 


88.0 

3.679 

3.482 

3.829 


90.0 

3.192 

3.040 

3.398 

2 

0.0 

2.368 

2.000 

1.777 


22.5 

2.323 

1.945 

1.693 


45.0 

2.212 

1.818 

1.559 


67.5 

2.121 

1.733 

1.635 


78.8 

2.087 

1.715 

1.713 


84.9 

2.111 

1.746 

1.783 


88.0 

2.142 

1.779 

1.830 


90.0 

1.958 

1.625 

1.646 
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Figure 4. Comparison of results for a surface crack emanating from the 
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Figure 9. Effect of r/t on the distribution of boundary-correction factors 
for a deep semi -el 1 iptical surface crack at the center of the 
notch root. 
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Figure 12. Effect of notch size on boundary-correction factors for 
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